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Abstract 

A search for selectron, smuon and stau pair production is performed with the 
data collected by the ALEPH detector at LEP at centre-of-mass energies up to 
209 GeV. The numbers of candidate events are consistent with the background 
predicted by the Standard Model. Final mass limits from ALEPH are reported. 
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1 Introduction 



The final results of searches for sleptons (£) with the data collected by the ALEPH detector 
at LEP at centre-of-mass energies up to 209 GeV are presented in this letter. These 
searches are interpreted in the theoretical framework of the minimal supersymmetric 
extension of the Standard Model (MSSM) with R-parity conservation and the 
assumption that the lightest neutralino is the lightest supersymmetric particle (LSP). 

Scalar leptons are pair produced at LEP through s-channel exchange of a Z or a 7 
giving rise to either £l£l or £r£r. Selectrons can also be produced through the t-channel 
exchange of a neutralino, resulting also in Ql^r production. Since the amount of mixing 
between the two slepton states (£l and £r) is proportional to the ratio of the lepton 
mass to the slepton mass, it is negligible for selectrons and smuons and may be sizeable 
for staus. As a consequence, the stau production cross section depends on the MSSM 
parameters A T , /1 and tan/3 via stau mixing, while the smuon production cross section 
depends only on the smuon mass. A dependence on M 2 , n and tan (3 is also present for the 
selectron cross section via the couplings with the neutralino in the t-channel production. 

Sleptons decay predominantly into their Standard Model (SM) partners and the 
lightest neutralino, £ ± — > P^Xi- The experimental signature is therefore a pair 
of oppositely-charged, same-flavour, acoplanar leptons (electrons, muons or taus), 
accompanied with missing energy carried away by the two undetected neutralinos. 

Results of previous slepton searches up to yfs = 202 GeV have been reported by 
ALEPH and by the other LEP collaborations ||. In this letter the analyses 

performed on the data taken in the year 2000 at centre-of-mass energies ranging from 
204 to 209 GeV are presented. The selections are mostly independent of the centre-of- 
mass energy except for an appropriate rescaling of the cuts with a/s when relevant. To 
derive the final results, all the data collected from 1997 to 2000 are used, with luminosities 
and centre-of-mass energies as shown in Table |1| 

Table 1: Integrated luminosities, centre-of-mass energy ranges and mean centre-of-mass energy 
values for data collected by the ALEPH detector from 1997 to 2000. 



Year 


Luminosity (pb x ) 


Energy range (GeV) 


<>/*> (GeV) 


2000 


9.4 


207 - 209 


208.0 




122.6 


206 - 207 


206.6 




75.3 


204 - 206 


205.2 


1999 


42.0 




201.6 




86.2 




199.5 




79.8 




195.5 




28.9 




191.6 


1998 


173.6 




188.6 


1997 


56.8 




182.7 
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This letter is organised as follows. In Section 0, a brief description of the ALEPH 
detector is given. The signal and background simulations are presented in Section [3[ In 
Section |, the selection criteria are described. Systematic uncertainties are discussed in 
Section [5]. The results are reported in Section |[ 

2 ALEPH detector 

The ALEPH detector is described in detail in Ref. ||. An account of its performance and 
a description of the standard analysis algorithms can be found in Ref. [[F]. 

In ALEPH, the trajectories of charged particles are measured with a silicon vertex 
detector, a cylindrical drift chamber, and a large time projection chamber (TPC). 
These detector components are immersed in a 1.5 T axial magnetic field provided by a 
superconducting solenoidal coil. Reconstructed charged particle trajectories are called 
good tracks if they are reconstructed with at least four space points in the TPC, a 
transverse momentum in excess of 200MeV/c, a polar angle with respect to the beam 
such that jcos^l < 0.95, and originating from within a cylinder of length 10 cm and 
radius 2 cm coaxial with the beam and centred at the nominal interaction point. 

The electromagnetic calorimeter (ECAL), placed between the TPC and the coil, is 
a highly-segmented, 22-radiation-length-thick sandwich of lead planes and proportional 
wire chambers. It consists of a barrel and two endcaps. It is used to identify electrons and 
photons by the characteristic longitudinal and transverse developments of the associated 
showers, supplemented for low momentum electrons by the measurement in the TPC of 
the specific energy loss by ionization. Photon conversions to e + e~ are identified as pairs 
of oppositely-charged electrons satisfying stringent conditions on their distance of closest 
approach and their invariant mass. The luminosity monitors (LCAL and SiCAL) extend 
the ECAL hermeticity down to 34 mrad from the beam axis. The hadron calorimeter 
(HCAL) consists of the iron return yoke of the magnet instrumented with streamer tubes. 
It provides a measurement of the hadronic energy and, together with the external muon 
chambers, efficient identification of muons by their characteristic penetration pattern. 

Global event quantities (such as total energy and missing momentum) are determined 
with an energy-flow algorithm which combines all the above measurements into charged 
particles (electrons, muons and charged hadrons), photons and neutral hadrons, called 
energy-flow particles in the following. Tau identification proceeds by clustering the energy- 
flow particles into two jets with the Durham algorithm ||. A jet is called a tau-jet 
candidate if it contains one or three good track(s) (or two if it contains an identified 
electron, to allow for asymmetric photon conversions) and if the jet invariant mass is 
smaller than 2GeV/c 2 . 
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3 Simulated samples 



The signal samples were generated using SUSYGEN , with final state radiation simulated 
with the PHOTOS JHJ package and tau decays simulated with TAUOLA [TT]. The slepton 
masses were varied in steps of 5GeV/c 2 (lGeV/c 2 ) for a mass difference AM with the 
lightest neutralino larger (smaller) than 10 GeV/c 2 ; the mass of the neutralino was varied 
steps of 5 GeV/c 2 . Event samples of all SM background processes relevant for the slepton 
search were also generated, i.e., dilepton production, 77 interactions with leptonic final 
states, W and Z pair production, and Zee and Wez/ production. The Bhabha process 
was simulated with BHWIDE |l2j, muon and tau pair production with KORALZ []13 1 , the 77 



processes with PH0T02 [III, WW production with KORALW |I3[ and the remaining four- 
fermion processes with PYTHIA fl6f . A detailed GEANT [17] simulation of the detector 
response was applied to both background and signal events. 



4 Description of the selections 



The final state topology depends on the slepton flavour. Selectrons and smuons are 
selected by requiring exactly two good tracks identified as electrons or muons and with 
opposite electric charge, while staus are selected by requiring two to six good tracks, 
clustered into two tau-jet candidates. The final state topology also depends on AM. 
Because the backgrounds are different for large values of AM (mostly W pair production) 
and small values of AM (mostly 77 interactions), the selections were optimized for two 
AM ranges. They are called "small-AM" and "large-AM" analyses in the following. The 
selections were optimized by minimizing the signal cross section expected to be excluded, 



on average, in absence of signal (the A95 prescription 0]). The selection criteria are 
summarized in Table for selectrons and smuons and in Table |3] for staus. The variables 
used in both the small- and large-AM analyses are defined below. 



• Acoplanarity <I> aC o 

The acoplanarity is defined as the angle between the two lepton momenta projected 
onto the plane transverse to the beam. An acoplanarity cut is used to reject events 
from e + e~ — > £ + £~ or 77 — > £ + £~ . 

• The p variable 

The lepton momenta are projected onto the transverse plane and the thrust axis is 
computed from the projected momenta. The p variable is defined as the scalar sum 
of the transverse momentum components of all energy-flow particles with respect 
to the thrust axis. A cut on p reduces the number of background events from the 
e + e~ — > t + t~ and 77 — > t + t~ processes. 

• Variables M V - 1S} p m [ SS , m i SS} m i SS , Prmiss an d PLmiss 

The total energy and momentum is constructed by adding energies and momenta 
of all energy- flow particles of the event, allowing a total visible mass M vis to be 
computed. The missing momentum p m i SS , identical in magnitude and opposite in 
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direction to the total momentum, defines the variables # m i SS and m i SS (polar and 
azimuthal angles of its direction), and pimiss andpLmiss (transverse and longitudinal 
components). Because the small- AM analysis is more sensitive to accidental 
double counting between low energy tracks and calorimeter clusters, two alternative 
methods are used to determine the total energy and momentum, (i) without the 
neutral hadrons (variables labelled "nH"); and (%%) without the neutral particles 
(variables labelled "tr"). The cuts on the visible mass and on the missing momentum 
(value and angle) reject background events from 77 processes. The cut on pTmiss is 
tightened when the missing momentum points to the vertical boundaries between 
the two halves of the LCAL. 

• Lepton momenta pi, p 2 and pji 

The variable pi (P2) is the momentum of the more (less) energetic identified electron 
or muon, if any. For massive neutralinos, upper cuts on M vis and pi are effective at 
rejecting leptonic WW events. 

• Activity at small angle E i2 and E 12 (H) 

Because 77 interactions often produce particles strongly boosted forward or 

Table 2: Selection criteria for the searches for acoplanar selectrons and smuons. The energies 
are expressed in GeV, the momenta in GeV /c and the masses in GeV /c 2 . Di-fermion final state 
events, four-fermion final state events and photon conversions are denoted by 2f, 4f and 7 conv., 
respectively. The sliding energy cut is explained in Section |6[ 



Selectron and smuon selection criteria 




Large AM 


Small AM 


Preselection 


Good tracks 


Two oppositely-charged, same-flavour, 
identified e or fi 


Anti-2f 


Acoplanarity 


$aco< 170° 


Anti-(2f+ 7 ) 


Neutral veto 


Cut applied 


Anti-(7 conv.) 


Acollinearity 


a > 2° 


Anti-77 


Energy 
within 12° 


£12 = 


E 12 (H) = 




Visible mass 


M vis >4 


M vis (nH) > 4 




P 


p > 2 


p>l 




Missing 


PTmiss > 3%a/s 


PTmiss(nH) > l%yfs 




momentum 


if | cos0 miss | < 0.26, 
then pxmiss > 5%a/s 


PTmiss(tl') > \%y/s 

| cos# miss (nH)| < 0.9 




Fisher variable 




Cut applied 




Momenta of the 


Pti,Pt2 > 0.5%v/i 


Anti-4f 


leptons 


pi < 46.5%v/i 


Pti < 10%v/i 




Visible mass 




M vis (nH) < 20% y/s 




Miss, momentum 




p miss (nH) < 10%v/i 




Sliding energy cut 
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backward, the requirement that no energy be detected at low angles is effective 
in suppressing this background. The variable E i2 gives the total energy measured 
within 12° of the beam axis in LCAL and SiCAL, and in the first cells (closest 
to the beam pipe) of the ECAL endcap calorimeters. The variable £12 (H) takes 
into account in addition the energy collected in the first cells of the HCAL endcap 
calorimeters. This cut introduces an inefficiency due to beam related background 
(Section ||). 

• Neutral veto 

Dilepton events with hard initial or final state radiation at large angle from incoming 
and outgoing leptons give rise to acoplanar lepton pairs, and to potentially large 
missing energy if the radiated photons go into poorly instrumented regions of the 
detector (e.g., ECAL module boundaries, overlap between barrel and endcaps). To 
reject these events, it is required that the invariant mass between any good track 
and any neutral energy-flow particle, outside a cone of half-angle 10° around the 
good track and with an energy in excess of 4 GeV, be smaller than 2 GeV/c 2 . 

Table 3: Selection criteria for the search for acoplanar staus. The energies are expressed in GeV, 
the momenta inGeV/c and the masses inGeV/c 2 . Di-fermion final state events, four-fermion 
final state events and photon conversions are denoted by 2f, 4f and 7 conv., respectively. The 
p T variable represents the smaller tau-jet momentum. The sliding energy cut is explained in 
Section |(| 



Stau selection criteria 




Large AM 


Small AM 


Pre- 


Good tracks 


Two to six tracks clustered in two r jet candidates 


selection 




p T x 206/v/i > 1.07 




Anti- 


Acollinearity 




a > 2° 


(7 conv.) 








Anti- 


Neutral veto 




yes 


(2f+7) 








Anti-2f 


Acoplanarity 


p x 206/v/i > ($ aco - 150°)/7.34 


Anti-77 


and p 








Energy 




E 12 = 




within 12° 








Visible mass 


M vis >4 




Missing 


Prmiss x 206/ y/s > 


PTmiss x 206/ v/i > 




momentum 


16.4 - M vis x 0.27 
PTmiss x 206/ v^i < 
(9 miss - 10°) x 1.24 


20.6-4.5 x p x 206/v/i 


Anti-4f 


Lept. momenta 


p x x 206/^ < 20.0 


pi x 206/v/i < 12.6 




Visible mass 


M vis < 80.0 






Sliding energy cut 
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Acollinearity a 

Single-photon events from e + e~ — > vv r y with a hard radiated photon, followed by 
a photon conversion into an e + e~ pair, are not rejected by the previous veto, but 
are characterized by a very small opening angle between the two tracks (called the 
acollinearity angle). An acollinearity cut is therefore effective at rejecting these 
events. 

Fisher variable 

In order to further reduce the large 77 background in the case of the small-AM 



selectron or smuon analyses, a Fisher discriminant [19| analysis is used. This method 
exploits the remaining modest difference among 77 events and the signal, taking 
into account the correlations between the variables M vis , pTmiss, PLmiss, P, Pi, P2, 
and a. 



5 Systematic effects 

The main systematic uncertainties on the background and signal expectations come from 
the statistics of the simulated samples (from 1.5 to 3.5%). A 2% systematic uncertainty 
was estimated for the lepton identification efficiency The inefficiency caused by the E12 
cuts is determined from events triggered at random beam crossings. The corrections 
for E 12 (E 12 (H)) used for the year 2000 data are -9.5% (-14.6%) at 205.2 GeV and 
—8.8% (—13.7%) at 206.6 and 208.0 GeV. The uncertainties on these corrections are less 
than 1%. 

The numbers of background events and the selection efficiencies are interpolated from 
the generated values to the mean energies of the year 2000 data. The contribution to 
the systematic uncertainty of this interpolation procedure is negligible. A comparison 
between data and expected backgrounds from SM processes is presented in Fig. [I] for the 
large-AM analysis before anti-four-fermion cuts and without lepton identification. At 
this level, agreement is observed between data and the expected backgrounds, dominated 
by W pair production. 



6 Results 

After the selection criteria described in Section |] are applied, the numbers of candidate 
events observed in the data and those expected from SM background processes are given 
in Table f| for the selectron and smuon searches and in Table [5] for the stau search. 

The search for sleptons is performed as a function of the hypothetical slepton and 
neutralino masses. A sliding energy cut is applied (not accounted for in Tables [| and |5|) 
which requires that the lepton momenta be in the range kinematically allowed for a 
signal with the assumed values of slepton and neutralino masses. For a two-body decay 
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Figure 1: Comparison between data (dots with error bars) and expected backgrounds 
(histograms) before anti-four-fermion and lepton identification cuts. For selectrons and smuons: 
(a) total visible mass (M v i s ) and (b) momentum of the leading lepton (pi). For staus: (c) total 
visible mass (M v ; s ) and (d) momentum of the leading lepton identified in the tau jet (pi). The 
accumulation at zero in (d) corresponds to events with no leptons in either tau jet. 
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Table 4: Numbers of candidate events observed in the year 2000 data (-/V b s ) and of background 
events expected from SM processes (iVbkg), in the searches for selectrons and smuons (before 
the sliding energy cut). The numbers in parentheses give the WW contribution. 



V* 

(GeV) 


e 


P> 


Small AM 


Large AM 


Small AM 


Large AM 




Nbkg 




Nbkg 


iVobs 


Nhkg 




Abk g 


205.2 





0.12 
(0.09) 


15 


14.0 
(12.2) 





0.12 
(0.07) 


8 


12.5 
(11.4) 


206.6 





0.19 
(0.16) 


24 


22.9 
(20.0) 





0.20 
(0.11) 


28 


20.6 
(18.8) 


208.0 





0.02 
(0.01) 





1.76 
(1.53) 





0.17 

(0.08) 


3 


1.59 
(1.44) 



Table 5: Numbers of candidate events observed in the year 2000 data (iV b s ) and of background 
events expected from SM processes (iVbkg), in the search for staus (before the sliding energy 
cut). The numbers in parentheses give the WW contribution. 



(GeV) 


r 


Small AM 


Large AM 


Small or Large AM 




Nhkg 


N ohs 


AW 




A^bkg 


205.2 


5 


6.1 
(5.0) 


1 


5.4 
(4.4) 


5 


7.6 

(6.3) 


206.6 


12 


10.1 

(8.3) 


10 


8.8 
(7.2) 


15 


12.5 
(10.3) 


208.0 


3 


0.80 
(0.68) 





0.70 
(0.57) 


3 


1.01 
(0.83) 



£ ± — > ^Xii the minimum and maximum values for the decay lepton energy are given by 




where is the beam energy and where the lepton mass is neglected. In the case of staus, 
the exact formula is used to account for the non-negligible mass of the tau. An event is 
compatible with the mass hypothesis (mi, m x o) if both decay leptons satisfy Eq. ([I]). Any 
given event may therefore be compatible with a wide mass range, especially for large- AM 
values. In the decay of a stau, the resulting tau lepton will itself give rise to at least one 
neutrino and the lower limit on the energy cannot be applied. Only the upper limit can 
be used, which increases even further the compatible region. For the selectron, smuon 
and stau searches, the numbers of candidate events and expected background events are 
displayed in Fig. ^| as a function of the slepton and the neutralino masses, after the sliding 
energy cut is applied. 
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(a) 






Figure 2: Numbers of selected events in the year 2000 data (left) and of expected background 
events from SM processes (right) in the search for (a) selectrons, (b) smuons and (c) staus. 
The density of shading reflects the number of events as indicated by the vertical scales. The 
kinematically accessible region at 209 GeV is indicated by the vertical dotted lines. The diagonal 
dotted lines show the boundary of the domain allowed in the MSSM for a neutralino LSP. 
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6.1 Cross section upper limits 



The data show no evidence for an excess of events compared to the estimated background. 
The dominant WW background theoretically and experimentally well under control, is 
subtracted to derive 95% C.L. upper limits on the slepton production cross section times 
the branching ratio into £xv This partial background subtraction is performed with the 



likelihood ratio method |20|. Data from previous years are combined in the likelihood 
to increase the sensitivity. The systematic uncertainties on the efficiencies are taken into 
account with the method of Ref. 1 12111 . 



The resulting bounds on the cross section rescaled at \/s = 208 GeV are shown in 
Fig. H Cross sections exceeding 0.1 (0.2) pb are excluded in a significant part of the plane 
(m^, m x o) for selectrons and smuons (staus). In the small- AM region, the limits are less 
stringent because the softer decay leptons cause a drop in the selection efficiency The 
limits are also less stringent for slepton masses close to the kinematic limit, where only 
part of the data contributes. 



6.2 Mass exclusion limits 

In the MSSM framework, the previous bounds allow limits to be set on the slepton masses 
as a function of the neutralino mass. To this end, cross sections and branching ratios are 
calculated with the program SUSYGEN ||. For selectrons and smuons, mixing is expected 
to be negligible and limits are derived under the conservative assumption that only £r£r 
production contributes. Limits are determined for staus in mixed and unmixed scenarios. 
Without mixing, conservative limits are set again under the assumption that only trTr 
production contributes. If the no-mixing assumption is relaxed, the most conservative 
limit on the mass of the lightest stau fi is obtained with a mixing angle Of which minimizes 
the production cross section, i.e., such that f\ decouples from the Z boson (Of ~ 52°). 

The unification of gaugino masses at the GUT scale is assumed for the computations 
of the slepton branching ratios and of the selectron production cross section. Branching 
ratios for the slepton decay are calculated for ji = — 200GeV/c 2 and tan/? = 2. The 
branching ratio into £ Xi is nearly 100% except for small neutralino masses, in which 
case the cascade decay into £ x% followed by x° ~~ > or Xil ma y become kinematically 
allowed. The limits are computed under the conservative assumption that the selection 
efficiency for decay channels other than £^ — > £ Xi is zero. 

The 95% C.L. expected and observed bounds on the masses of selectrons, smuons and 
staus as a function of the neutralino mass are displayed in Fig. [|, together with the effect 
of cascade decays. For staus, the limit obtained with mixing angle Of ~ 52° is also given. 
The observed stau mass limit is less stringent than expected because a small excess of 
events was observed in 1999 ||. 

The 95% C.L. bounds on the slepton masses for AM > 15GeV/c 2 are summarized 
in Table |6], with the hypothesis BR(£ ± — > £^Xi) = 1- For smuons and staus, these limits 
are independent of the MSSM parameters. For selectrons, complementary analyses and a 
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(a) 





Figure 3: Upper limit at 95% C.L. on the observed (left) and expected (right) production cross 
section at yfs = 208 GeV, times the branching ratio into £xi f° r ( a ) selectrons, (b) smuons and 
(c) staus. 
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Figure 4: Excluded regions at 95% CL. in the versus m^o plane from slepton searches. The 

observed (shaded area) and expected (solid curve) limits are given for BR(^ =t — > € X?) = 1- For 
selectrons, /i = — 200GeV/c 2 and tan/3 = 2 are assumed. The dashed curves give the observed 
limits when slepton cascade decays are taken into account. For staus, the dotted curve gives the 
most conservative limit, obtained for Of ~ 52° and with cascade decay taken into account. For 
the dashed and dotted curves, BR (^ ± — > ^x?) is computed with the values \i = — 200GeV/c 2 
and tan (5 = 2, and a zero efficiency for the selection of the topologies arising from cascade 
decay is assumed. The kinematicaily accessible region at 209 GeV is indicated by the vertical 
dash-dotted lines. The diagonal dash-dotted lines show the boundary of the domain allowed in 
the MSSM for a neutralino LSP. 
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Table 6: Lower limits at 95% C.L. on the slepton masses for AM > 15GeV/c 2 , with the 
hypothesis BR^ 1 * 1 — > I Xi) = 1- For staus, the bound corresponding to the minimum cross 
section (as explained in the text) is also given. In the case of selectrons, the limits are given for 
fj. = -200GeV/c 2 and tan/3 = 2. 



Channel 


m(£) > 
obtained 


m{l) > 
expected 


P>R 
TR 

f(min) 


95 GeV/c 2 
88 GeV/c 2 
79 GeV/c 2 
76 GeV/c 2 


96 GeV/c 2 
87 GeV/c 2 
83 GeV/c 2 
81 GeV/c 2 



scan on the relevant MSSM parameters allow an absolute mass lower limit to be extracted. 



The result of such a study is presented in Ref. [22 . 



7 Conclusions 

Searches for scalar lepton pair production have been performed in the data sample 
collected by the ALEPH experiment at LEP2 at centre-of-mass energies up to 209 GeV. 
The numbers of candidate events observed are consistent with the background expected 
from Standard Model processes. 

In the framework of the MSSM, 95% C.L. mass exclusion regions have been obtained 
for selectrons, smuons and staus in the plane (ra^, m x o). In particular, mass lower limits 
for smuons and staus are set at 88 and 76GeV/c 2 , respectively, for a mass difference 
between the slepton and the lightest neutralino in excess of 15GeV/c 2 and for a slepton 
decay branching ratio into l\\ of 100%. These limits are independent of the MSSM 
parameters. Under the same assumptions, and with fi = — 200GeV/c 2 and tan/3 = 2, a 
lower limit of 95GeV/c 2 is set on the selectron mass. 
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